Optical Flux Lattices for Ultracold Atomic Gases 
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We show that simple laser configurations can give rise to "optical flux lattices" , in which optically 
dressed atoms experience a periodic effective magnetic flux with high mean density. These potentials 
lead to narrow energy bands with non-zero Chern number. Optical flux lattices will greatly facilitate 
the achievement of the quantum Hall regime for ultracold atomic gases. 



PACS numbers: 67.85.Hj, 67.85.-d, 37.10.Vz 

One of the most important techniques in the ultra- 
cold atom toolbox is the optical lattice[l|: a periodic 
scalar potential formed from standing waves of light. 
Optical lattices are central to the use of atomic gases 
as quantum simulators, and allow the exploration of 
strong-correlation phenomena related to condensed mat- 
ter systems [1, Their usefulness derives from the fact 
that the scalar potential has a lengthscale, set by the 
optical wavelength A, that is similar to the typical inter- 
atomic spacing. 

Largely separate have been the experimental develop- 
ments of effective vector potentials, representing the cou- 
pling of a charged particle to a magnetic field. A uniform 
magnetic field can be simulated for neutral atoms using 
rotationjl,, 4, 5^. Gauge fields leading to an effective mag- 
netic field can also be generated by optical dressing [6[. 
These techniques have allowed experimental studies of 
quantized vortices in condensed Bose and (paired) Fermi 
gases 111 0, Q-Q- interesting regime of strong corre- 
lation, with connections to the fractional quantum Hall 
effect, is expected when the effective magnetic flux den- 
sity is sufficiently high that the lengthscale n^"^/^ is 
comparable to the interatomic spacing|5i [lo| . However, 
the magnetic flux densities achievable using existing tech- 
niques are relatively small for large systems d, [ll|, so 
this strongly correlated regime occurs at very low densi- 
ties when the gas is weakly interacting and highly sus- 
ceptible to uncontrolled perturbations [Sl. 1 11 1 . 

In this paper we describe simple laser configurations 
that use optical dressing to generate what we refer to 
as "optical flux lattices". While a conventional optical 
lattice imprints a periodic scalar potential, an optical 
flux lattice also imprints a periodic magnetic flux density 
with non-zero mean, and large magnitude fi^ ~ 
We emphasize that the magnetic flux density is a con- 
tinuous function of position, so these potentials are dis- 
tinct from proposed ways to imprint gauge potentials 
onto deep optical lattices which apply only in the tight- 
binding limit[12Hl4l|: furthermore, optical flux lattices re- 
quire only a small number of lasers, so are much easier 
to implement than these tight-binding proposals. We 
show that optical flux lattices lead to narrow bands with 
non-zero Chern number [l5j|. In particular the lowest en- 
ergy band is topologically equivalent to the lowest Lan- 
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dau level. Since the lengthscale ^ A is similar to 

the typical interatomic spacing, optical flux lattices will 
allow the study of quantum Hall physics at high densities 
where interaction energy scales are large. 

We consider an atom moving in optical fields within 
the rotating wave approximation, with Hamiltonian 
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where V is the energy scale of the optical potential, of 
dimensionless form M{r). We focus on two-level systems 
and write 



M ^ M{r) ■ a 



+ iMy 
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with a the vector of Pauli matrices. Additional scalar po- 
tentials can be added by conventional optical lattices; for 
simplicity, we neglect these here. The off-diagonal terms 
[Mx^y) arise from the Raman coupling that effects the 
interspecies conversion. The diagonal term {Mz) repre- 
sents a species-dependent potential. One possible imple- 
mentation of the two- level system is with the groundstate 
and long-lived excited state of an alkaline earth atom or 
ytterbium 1J| , in which case can be generated by a 
laser at an "anti- magic" wavelength, Aam- As we shall 
describe below, optical flux lattices can be formed by a 
standing wave at Aam and three travelling waves of the 
Raman laser [3l . We shall restrict attention to quasi two- 
dimensional (2D) systems with r = {x,y). The resulting 
flux lattices are readily adapted to 3D, with net flux along 
one direction. 

The emergence of an effective gauge potential is best 
understood when the kinetic energy is small compared to 
the energy spacing of the (local) dressed states, obtained 
from the eigenvalues of VM. The atom then moves 
through space adiabatically, staying in a given dressed 
state. The adiabatic limit is always valid for sufficiently 
large V provided the spectrum of M is non-degenerate. 
Assuming this to be true (as shall be verified below for 
the cases of interest), we consider the adiabatic motion 
in a normalized dressed state 
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Projecting ([T]) onto adiabatic motion on the state ^ 
leads to an effective Hamiltonian with both a scalar and 
a vector potential, the latter given byQ 

qA^ ih{^\V^) (4) 

for effective charge q. The number density of magnetic 
flux quanta perpendicular to the xy-plane is therefore 

For optical fields of wavelength A, it is natural to as- 
sume that the vector potential ^ is smoothly varying 
with \qA\ ^ h/ X. Then, the maximum mean flux density 
in a region of space of sides L^, Ly 3> A may be found 
by applying Stokes' theorem: J n^dPr = h^LxLy — 
[q/h) § A - dr < {Lx + Ly)/X, leading to < 1/LX with 
L — min(i^, Ly\ All existing proposals for optically in- 
duced gauge fields in the continuum [ij, and the scheme 
implemented in are of this form with the scale L 

set by the width of the cloud. Since, typically, i » A, 
this leads to relatively small flux density, ~ 1/LX. 

Although apparently very general, these considerations 
neglect the fact that smoothly varying optical fields can 
induce singularities in qA. These singularities depend 
on the gauge used for ([3]), and cause no singularities in 
gauge-invariant properties. Such issues arise whenever 
a U(l) gauge field has non-zero flux through a closed 
manifold, notably leading to Dirac strings for a magnetic 
monopole. For the optical flux lattices we propose here, 
there is a net flux through a unit cell which (due to the 
spatial periodicity) has the topology of a torus. We avoid 
technical difhculties of the gauge-dependent singularities 
by defining the local Bloch vector 

n{r) = ($(r)|^|$(r)) (6) 

for which n ■ ft — 1 . The flux density is then 

^ --^e^jkef^^n^dfj^njdi^nk ■ (7) 

OTT 

This is (minus) the "topological density" of the map from 
position space, r = {x,y), to the surface of the Bloch 
sphere, n[l7[. The number of fiux quanta through a re- 
gion A is J^ncj,(Pr = U/Att where is the solid angle 
that region A maps to on the Bloch sphere. Thus, each 
time the Bloch vector wraps the surface of the sphere 
corresponds to one magnetic flux quantum. Optical fiux 
lattices are spatially periodic configurations for which the 
Bloch vector wraps the sphere an integer number, iV^, 
times in each unit cell. The lattice vectors ai and 02 are 
both of order the optical wavelength A, so the mean flux 
density is of order ^ N^/X^ which is large. We focus 
on two cases of high symmetry which achieve this goal. 
Square Lattice: Consider the optical coupling 

Mgq = cos{K,x)ax -f cos(Ky)iTj, -I- sin(Ka;) sm{K,y)az (8) 




FIG. 1: Properties of the lower energy dressed state for the 
square optical flux lattice ((8|, as a function of position in the 
unit cell, (a) Local Bloch vector, as represented by the vector 
(n^, Uy) and contours of (light shading for < 1 and dark 
for Uz > 1). (b) The local flux density is maximal on the 
lines X — ±a/4 and y — ±a/4 and vanishes at four points. 
The lattice sites and nearest-neighbour hopping in the tight- 
binding limit are indicated by dark circles and dashed lines. 

where k = lixja. This has square symmetry with 
a\ = (a, 0),a2 = (0,a). Achieving this high symme- 
try in experiment may involve tilting the lasers out of 
the xy-plane to tune the periods of the Raman and the 
species-dependent fields [l8|. The eigenvalues are non- 
degenerate at all positions, so the dressed states admit 
an adiabatic limit. Figllja) shows that, for the lower en- 
ergy dressed state, the (n^^ny) components of the Bloch 
vector exhibit two vortices and two anti-vortices in the 
unit cell. These vortices lead to gauge-dependent singu- 
larities in the vector potential. However, the Bloch vec- 
tor varies smoothly, with = ±1 at the vortex cores, 
in such a way that all four of these regions wrap the 
sphere in the same sense and contribute a fiux density 
of the same sign. The fiux density is shown in Figlljc). 
It is everywhere non-negative and has total fiux = 2 
in the unit cell. (This may be seen by noting that the 
two vortices and two antivortices cause n to wrap the 
Bloch sphere twice.) The flux density is not constant, 
and vanishes at four points in the unit cell. These four 
points coincide with the locations at which the adiabatic 
energy is minimum, FigHIb), so form the lattice sites in 
the tight- binding limit. 

Triangular lattice: An optical flux lattice with trian- 
gular symmetry is generated by 

Mtri = cos(r • Ki)(Jx + cos(r • K2)a-y + cos(r • K3)(Jz (9) 

where Ki = (1,0)k,K2 = (l'"^)'^ «3 = «i — 

K2, with K = 47r/(\/3a), giving lattice vectors ai = 
(V3/2,-l/2)a and = (0, l)a. A gain, the eigenval- 
ues are non-degenerate at all positions. The properties 
of the lower energy dressed state are shown in Figl2] The 
fiux density is maximum on the sites of a honeycomb lat- 
tice, and vanishes on the triangular lattice dual to this. 
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FIG. 2: (a) Bloch vector and (b) flux density for the lower 
energy dressed state of the triangular optical flux lattice ((8|. 
The local minima of the adiabatic energy are at the points 
where = 0, forming a triangular lattice in the tight-binding 
limit (dark circles and dashed lines). 



The tight-binding limit involves tunneling between the 
sites of this dual triangular lattice, FiglTJb). 

The above optical potentials (I8I9P can be readily gen- 
erahzed to many other cases with non-zero mean flux. 
(There are also many cases with zero mean, but non-zero 
local flux density.) The central requirements to generate 
an optical flux lattice are threefold. Firstly, the Raman 
laser (M^.y) must generate optical vortices. A 2D lat- 
tice of optical vortices can be formed from a minimum of 
three travelling waves The resulting optical field is 
periodic, with an equal number of (single-winding) vor- 
tices Nv and antivortices A'^av in a unit cell 20]. Secondly, 
the species-dependent potential (M^) must be non-zero 
at the cores of these vortices, such that there is no de- 
generacy of the dressed states at these points. A small 
non-zero causes the cores of the vortices to have 
the topology of "merons" [l7j . in which n{r) sweeps over 
half of the Bloch sphere. For a given meron, the sign 
of Mz at its core times the sign of its vorticity deter- 
mines whether it contributes -1-1/2 or —1/2 a flux quan- 
tum. The total number of flux quanta through the unit 
cell is Nd, — N.f — N^., where A^^, is the number of 

V V ? V / av 

vortices/antivortices at which is positive. Thus, the 
third requirement for a nonzero mean flux is that 
varies in space such that A'^ ^ A^^!^. 

We have explored the properties of optical potentials 
generated by simple laser patterns. An optical flux lat- 
tice can be generated using just five travelling waves: 
three travelling waves of the Raman laser (Mx^y) to ef- 
fect the vortex lattice, and a standing wave of the species- 
dependent potential (M^). [One such example is to re- 
move one of the four travelling waves from the Raman 
coupling in In all cases the local fiux density is inho- 



mogeneous in space, in some even changing sign. Indeed, 
one can show that, for smoothly varying optical fields, the 
flux density must have at least A'v -t- A^av zeroes in the 
unit cell[2l| The above cases (j8l9|) have non-negative fiux 
density with the minimum number of zeroes. For three-, 
or more-, level systems, an optical flux lattice can have a 
flux density that nowhere vanishes. We have examples of 
optical potentials that lead to such cases. However, these 
require more involved laser configurations, so we do not 
pursue this here. 

Having determined the properties of the optical fiux 
lattices in the adiabatic limit, we now turn to describe 
their bandstructures, obtained from the eigenvalues of 
dl]). The laser potentials M^q ^ and Mtri © are clearly 
invariant under translations by the respective lattice vec- 
tors ai_2- In fact, they enjoy higher translational sym- 
metry, being invariant under the unitary transformations 



To 



(10) 



which effect translations by ^ai,2 ^md rotations in spin- 
space. These operators do not commute, but satisfy 



T2T1 = -T1T2 



(11) 



This indicates that they represent magnetic translations 
around a region of space (enclosed by |ai and ^0,2) that 
contains 1/2 a flux quantum. As is conventional in sys- 
tems with magnetic translation symmetry 15| , we define 
a magnetic unit cell that encloses an integer number of 
fiux: we choose ai, 02/2. Writing the eigenvalues of the 
associated (commuting) translation operators and T2 
as e'^'' "'^ and e*'"'"'^/^ defines the Bloch wavevector k and 
the associated Brillouin zone. The additional symmetry 
Ti and the condition (ITU cause the energy spectrum Ek 
for all bands to be invariant under fc-a2/2 — )■ fe-a2/2±7r. 

For the square optical fiux lattice ([8]), a solution of 
the bandstructure shows that the lowest energy band 
does not overlap any higher band for V ^ O.lfi^K^/m. 
The Chern number [l5| of this band is 1, the sign being 
reversed under an odd number of sign changes to the 
terms in (|8]). Thus, the lowest energy band is topologi- 
cally equivalent to the lowest Landau level of a charged 
particle in a uniform magnetic field. It is instructive to 
consider the bandstructure for V ^ Tt?k^ /m, when the 
variation in the adiabatic energy is dominant, and the 
low energy bands are well described by a tight-binding 
modelQ. The minima of the adiabatic potential form 
a square lattice, Figdjb). Nearest- neighbour hopping on 
this square lattice leads to a model in which each plaque- 
tte encloses 1/2 a fiux quantum. The magnetic unit cell 
contains two lattice sites, so there are two ti ght- binding 
bands. The bands touch at two Dirac points[2^, so one 
can speak only of the Chern number of the two bands to- 
gether. This total Chern number is zero, consistent with 
the fact that this nearest-neighbour tight-binding model 
is time- reversal symmetric 23[ . In the physical model. 
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FIG. 3: Lowest energy bands for the triangular flux lattice 
^ in the nearest-neighbour tight-binding limit (for uniform 
spacing of — 27r/\/3 < k^a < 27r/%/3)- The energy Etb is 
relative to the atomic limit, in units of the nearest-neighbour 
hopping. The bands have Chern numbers ±1. 

with mVlfP'K^ large but finite, time-reversal symmetry 
is broken by next nearest-neighbour hopping across di- 
agonals of the square lattice. This leads to closed loops 
around plaquettes which contain 1/4 of a flux quantum. 
This perturbation acts to split the bands at the two Dirac 
points, and the two bands acquire Chern numbers of ±1. 

The bandstructure of the triangular optical flux lat- 
tice ^ has the same qualitative properties, the lowest 
energy band having a Chern number of 1. In this case, 
time-reversal symmetry is broken even in the tight bind- 
ing limit with nearest-neighbour hopping. The energy 
minima are at the sites of a triangular lattice, Figlljb), 
the elementary plaquettes of which enclose 1/4 of a flux 
quantum. The energy spectrum of the resulting tight- 
binding model, shown for a convenient gauge in Fig|31 
has two narrow bands that are well separated in energy 
and have Chern numbers of ±1. 

Optical flux lattices will allow experiments on ultra- 
cold gases to explore many very interesting phenomena. 
Since they lead to a lowest energy band with non-zero 
Chern number, non-interacting fermions filling this band 
(with one fermion per magnetic unit cell) will exhibit 
the integer quantum Hall effect. Signatures of the re- 
sulting chiral edge state could be observed in the density 
excitations (collective modes), which will rotate with a 
handedness determined by the sign of the Chern number. 
The square lattice in the nearest-neighbour tight binding 
limit also offers the possibility to study fermionic Dirac 
physics. Within mean-field theory, interacting bosons 
loaded into the chiral band will develop vortex lattices 
with very high fiux density. These typically break trans- 
lational symmetries of the lattice [24l|. Owing to the very 
high flux density, it should be possible to reach a regime 
where the 2D boson density is comparable to the mean 
flux density, ~ 1/A^, where strongly correlated frac- 
tional quantum Hall states of bosons [si [l^l or related 
states on lattices[2^ can appear. A leading candidate 



is the V = \/2 bosonic Laughlin state on the triangu- 
lar lattice ([9]), for which the lowest energy chiral band is 
narrow and well-separated from higher bands. There is 
1/2 a flux quantum per lattice site, so the Laughlin state 
appears at 1/4 filling. It will be interesting also to ex- 
plore strong correlation phenomena in 3D settings, with 
an optical flux lattice providing net flux in one direction. 
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